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I.   INTRODUCTION 

Setting  production  targets  for  geographically  dispersed  production  sites 
is  a  common  problem  in  large  manufacturing  companies.    Although  it  is  a 
routine  decision,  it  encompasses  important  but  sometimes  subtle  tradeoffs 
that  have  a  direct  impact  on  corporate  performance.   The  location  of 
production  clearly  affects  the  cost  of  production,  because  the  cost  of  inputs,  as 
well  as  taxes  and  other  financial  factors,  may  show  wide  geographic  variation. 
It  also  affects  the  cost  of  distributing  products  to  customers,  which  is  largely  a 
function  of  distance.   Service  levels  are  also  strongly  related  to  the  proximity 
of  the  production  point  to  customers,  and  to  such  factors  as  the  capacity  of  the 
plant. 

Production  targets  also  set  the  stage  for  lower  level  manufacturing 
decisions,  such  as  production  scheduling  and  inventory  management.    More 
generally,  targets  largely  determine  the  level  of  utilization  for  manufacturing 
sites.   At  the  most  extreme,  the  allocation  of  production  may  call  for  a  site  to 
be  shut  down  or  mothballed. 

This  paper  describes  a  production  allocation  modeling  system 
(PRAMSYS)  that  we  have  done  for  the  ArON  corporation,  a  major  producer 
of  industrial  gases  (oxygen,  nitrogen,  argon),  with  numerous  manufacturing 
sites  and  customers  throughout  the  United  States.   The  system  has  been  in 
active  use  as  a  routine  planning  tool  for  more  than  a  year  in  the  company's 
Mid-Atlantic  region,  and  installations  in  other  regions  are  underway. 

PRAMSYS  is  a  large  Mixed  Integer  Programming  (MIP)  model 
embedded  in  a  menu-driven  interface  for  data  preparation  and  report 
generation.    The  immediate  purpose  of  PRAMSYS  is  to  minimize  combined 
regional  manufacturing  and  distribution  costs  over  a  planning  period  of  one 
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to  three  months.    More  generally,  corporate  planners  use  it  each  month  to 
allocate  individual  customer  demands  to  geographically  dispersed  production 
sites.   PRAMSYS  also  optimally  allocates  idle  time  to  sites  in  keeping  with 
complex  relationships  between  production  costs  and  capabilities. 

PRAMSYS  is  of  general  interest  beyond  the  immediate  application 
because  it  demonstrates  the  value  of  optimization  in  a  mature  industry 
where  conventional  wisdom  might  lead  one  to  expect  opportunities  for  cost 
reduction  to  be  limited.   In  particular,  it  demonstrates  how  an  integrating 
model  can  be  used  to  bring  the  company's  technical  expertise  in  production 
and  process  engineering  to  bear  on  larger  strategic  issues. 

The  PRAMSYS  project  also  illustrates  how  a  model  and  an  application 
evolve  together  over  the  course  of  a  project  through  an  interplay  between 
practical,  computational,  and  theoretical  considerations.    In  this  case  the 
model  became  both  more  correct  and  simpler  -  a  happy  but  perhaps 
fortuitous  outcome  which  is  by  no  means  the  rule  with  complex  modeling 
applications.    Among  other  things,  this  evolution  illustrates  the  subtle 
difficulties  that  can  arise  when  practitioners  focus  too  closely  on  their 
mathematical  abstractions  and  lose  sight  of  the  practical  reality  behind  their 
models. 

II.    INDUSTRY  BACKGROUND 

Production  of  industrial  gases  is  in  many  ways  the  quintessential 
mature  manufacturing  industry.   The  process  of  cryogenic  distillation  by 
which  air  is  separated  into  gaseous  and  liquid  elemental  fractions  has  been 
known  for  over  eighty  years.   Competing  producers  now  operate  capital 
intensive  plants  with  similar  intrinsic  thermodynamic  efficiencies;  no  radical 


breakthroughs  in  production  technology  are  to  be  expected.   Air,  the  sole  raw 
material,  is  free  and  does  not  vary  appreciably  in  quality.   Nor  is  there  much 
scope  for  product  differentiation  -  except  for  special  applications  where 
extreme  purity  is  essential,  all  liquid  nitrogen  is  very  much  the  same. 

Despite  this  stability  on  the  supply  side,  however,  the  markets  for 
industrial  gases  are  changing,  largely  in  response  to  structural  changes  in  the 
national  and  world  economy.   Demand  for  liquid  and  gaseous  oxygen  was  for 
many  years  the  driving  force  of  the  industry.   In  recent  years  this  demand  has 
been  declining  as  the  centers  of  basic  industries  such  as  steelmaking  shift 
offshore.   On  the  other  hand,  demand  for  liquid  nitrogen  is  increasing  for  use 
in  food  preparation,  enhanced  oil  recovery,  and  other  domains  where  a 
combination  of  very  low  temperature  and  chemical  inertness  is  essential. 
These  changes  have  also  led  to  a  shift  in  the  location  of  demand,  away  from 
the  older  midwestern  industrial  centers. 

The  result  has  been  to  alter  the  prevailing  premises  and  operating 
procedures  in  the  industry.   The  company  is  no  longer  principally  an  adjunct 
of  stable  larger  industries  and  cannot  afford  to  operate  as  if  it  were.  This  shift 
in  the  conditions  underlying  competition  in  the  industry  raises  hazards 
where  for  decades  there  had  been  stability.   It  also  opens  up  new  opportunities 
for  a  company  that  can  adapt  to  the  new  conditions. 

III.   PROJECT  BACKGROUND 

The  PRAMSYS  project  originated  in  a  general  desire  on  the  part  of 
upper  management  at  ArON  to  bolster  competitive  position  through  better 
operation  of  the  company's  production  and  distribution  system.    Delivered 
cost  is  one  of  the  primary  determinants  of  competitive  advantage  in  this 


industry  (the  other  being  customer  service.)  The  two  primary  elements  of  cost 
that  are  subject  to  control  over  the  short  and  medium  term  are  distribution 
and  production.   Production  costs  are  generally  larger,  but  distribution  costs 
are  still  significant  -  on  the  order  of  35%.   It  was  therefore  natural  that 
attention  should  have  focused  at  first  on  reducing  each  of  these  costs 
independently  of  the  other,  particularly  since  such  an  effort  meshed  with  the 
prevailing  division  of  responsibilities  under  the  company's  Distribution  and 
Production  functions.   The  use  of  formal  models  closely  paralleled  this 
division. 

For  the  purposes  of  production  and  distribution  planning,  ArON 
groups  its  customers  and  production  sites  into  several  large  regions.    In 
principle,  any  site  in  a  region  can  serve  any  customer  in  that  region,  provided 
it  makes  the  product  demanded  by  the  customer.    Within  a  region,  known  or 
predicted  customer  demands  are  assigned  to  a  site  through  a  monthly 
planning  cycle.   These  demands  can  then  be  aggregated  into  production 
targets  for  each  product  at  each  site. 

In  practice,  the  regional  distribution  department  assigned  customers  to 
sites,  since  it  was  they  who  managed  the  physical  shipment  of  product  to 
customers.    In  making  this  decision,  the  department  made  heavy  use  of  a 
model  based  on  network  algorithms  which  had  been  developed  to  minimize 
distribution  costs.   This  naturally  tended  to  favor  assignment  of  customers  to 
the  nearest  site,  with  out  taking  fully  into  account  the  cost  of  production  at 
the  site.   The  production  process  and  its  economics  were  simply  too  complex 
to  be  represented  well  in  such  a  model.   In  a  very  real  sense,  the  model  served 
only  to  formalize  the  standard  operating  procedures  of  the  distribution 


department  and  to  optimize  the  allocation  of  production  with  reference  to 
interests  and  capabilities  of  that  department. 

The  production  function  also  had  in  place  a  very  successful  program  to 
improve  the  efficiency  of  production  sites.    A  major  element  of  this  program 
was  a  set  of  data  gathering  procedures  and  modeling  software  called  the  Site 
Process  Optimization  Protocol  (SIPOP).   ArON  chemical  and  process 
engineers  had  developed  SIPOP  for  use  at  each  site  to  determine  how  the  site 
should  be  operated  to  minimize  electric  power  demand.   Given  a  set  of  target 
production  rates  for  each  product,  SIPOP  calculates  the  minimum  power 
demand  to  produce  at  those  rates,  and  indicates  how  the  plant  should  be 
configured  and  operated  to  attain  this  minimum.    Although  SIPOP  performs 
this  localized  optimization  rapidly  and  very  accurately,  it  could  not  in  itself 
determine  what  the  target  production  rates  should  be,  since  these  directly 
reflect  higher  level  decisions  about  the  allocation  of  customer  demand  to  the 
site. 

As  we  reviewed  ArON's  procedures  and  planning  tools  it  became  clear 
that  cost  reductions  in  production  and  distribution  would  be  at  best  be 
haphazard,  if  not  illusory,  unless  they  were  achieved  in  concert.   Conspicuous 
by  its  absence  was  the  ability  to  plan  both  production  and  distribution 
activities  within  a  single,  comprehensive  framework  to  achieve  the  greatest 
overall  cost  reductions. 

This  kind  of  coordinated  planning  looked  to  be  relatively  untapped 
area  in  which  the  company  could  distinguish  itself  from  its  competitors.   In  a 
competitive  industry  such  as  industrial  gases  cost  reductions  of  even  one  or 
two  percent  can  be  extremely  important,  translating  into  much  larger 
increases  in  profit. 


IV.   MODEL  DEFINITION 

Air  separation  sites  produce  gaseous  and  liquid  air  fractions.  Gases  are 
distributed  by  pipeline  to  customers  located  near  the  site.   Liquids  (oxygen, 
nitrogen,  argon)  are  distributed  by  truck  or  railroad  tanker.   There  are  no  joint 
deliveries;  in  fact,  each  vehicle  is  dedicated  to  a  single  product.   This  means 
that  the  distribution  system  -  and  costs  -  for  each  product  are  linked  only 
indirectly,  through  joint  production  at  the  sites. 

The  distribution  component  of  PRAMSYS  is  thus  represented  as  a  set 
of  simple,  single  commodity  arc  networks  (one  for  each  product)  linking 
production  (or  external  supply)  points  to  customers.   In  general,  any  site  can 
deliver  to  any  customer,  but  if  shipment  from  a  given  site  to  a  given 
customer  is  undesirable  or  impossible,  then  the  corresponding  arc  may  be 
omitted  from  the  network.   Unit  transportation  cost  between  a  site  and  a 
customer  reflects  the  distance  between  the  two  points,  and  perhaps  the 
intervening  geography.   The  costs  used  in  PRAMSYS  are  derived  from 
historical  data  and  were  already  in  use  for  distribution  planning. 

Most  of  the  structure  of  the  model  emerges  from  the  representation  of 
the  production  sites.   The  complexity  of  this  representation  stems  from  a 
variety  of  related  factors,  among  them  joint  production,  electricity  contracts, 
and  shut-down  operation. 

A.  joint  Production 

A  site  produces  products  jointly  from  the  same  production  process  -  up 
to  five  products  at  once.   A  product  can  be  produced  at  any  rate,  within  upper 


and  lower  limits  that  depend  upon  the  site,  how  the  site  is  configured,  the 
product,  and  the  rates  at  which  other  products  are  being  produced. 
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Figure  1 


The  (instantaneous)  power  demand  of  the  site  is  an  increasing  function 
of  production  rates  for  all  products,  with  strong  cross  terms,  particularly  for 
liquid  products.   There  are  compelling  theoretical  and  empirical  reasons  to 
believe  that  the  surface  is  convex,  but  no  closed  form  for  the  function 
KW=f  (Pl,...Pn)   is  known. 
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Figure  2 
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Data  that  defines  this  surface  in  PRAMSYS  is  ultimately  derived  from  SIPOP, 
although  SIPOP  was  not  developed  for  this  purpose.   The  methodology  of 
SIPOP  is  in  fact  radically  different  from  that  of  PRAMSYS.  SIPOP  uses 
random  search  methods  to  find  the  minimum  power  point.    This  is  standard 
practice  in  chemical  engineering  and  process  modeling,  where  the  complexity 
and  non-linearity  of  the  underlying  processes  make  gradient  search  methods 
difficult  to  implement  and  cumbersome  to  use  (see  Martin  (1982),  Wang 
(1978)). 

B.  Electricity  Contracts 

Virtually  the  only  variable  production  cost  is  the  cost  of  electricity  used 
to  run  compressors  and  liquefiers,  so  that  production  cost  is  very  closely  tied 
to  the  site's  use  of  electricity.   Because  KW  demand  is  a  function  of  the 
production  rates  of  all  products,  a  decision  to  assign  a  customer's  demand  for, 
say,  LOX  (liquid  oxygen)  to  site  A  therefore  alters  the  cost  of  both  LOX:  and  of 
LN  (liquid  nitrogen)  at  that  site,  even  though  the  production  rate  for  LN 
remains  the  same. 

But  production  cost  at  a  site  is  not  strictly  a  matter  of  thermodynamic 
efficiency.   Sites  are  such  major  consumers  of  electricity  that  energy  and 
power  costs  are  governed  by  special  contractual  terms  that  are  often  quite 
complex  and  that  differ,  sometimes  radically,  from  site  to  site. 

One  typical  contractual  feature  is  that  the  site  is  charged  both  for  energy 
(KWH)  consumption  and  for  maximum  (instantaneous)  power  (KW) 
demand  during  some  contract  billing  period.   These  costs  are  roughly  of  the 
same  magnitude,  although  energy  costs  tend  to  be  higher. 
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Under  most  contracts  the  unit  cost  of  energy  varies  discontinuously  by 
time  of  day.   Figure  3  depicts  a  situation  in  which  the  day  is  divided  into  on- 
peak,  mid-peak,  and  off-peak  hours.   Energy  charges  are  highest  during  the 
on-peak  hours,  lowest  during  off-peak,  and  take  on  an  intermediate  value 
during  midpeak.   The  relative  proportion  of  on-,  off-,  and  mid-peak  periods 
in  a  weekday,  weekend  day,  and  holiday  may  all  be  different.  Any  period  type 
may  be  absent  from  any  day  type. 


$/KWH 


on 

mid 

mid 

off 

off 

0  —  —  —  -  hour  in  day  —   —  —   —  24 
Figure  3 
Time  of  Day  Energy  Pricing 


The  charge  on  power  demand  ($/KW)  may  also  be  different  in  on,  off, 
and  mid-peak  periods,  although  not  necessarily  in  proportion  to  energy 
charges.  Often  demand  charges  are  only  incurred  during  certain  periods,  e.g., 
on-peak. 

Under  such  contractual  terms  there  is  a  strong  incentive  to  produce  at 
higher  rates  during  off-peak  periods,  when  energy  and  power  are  cheaper,  and 
to  throttle  back  during  more  expensive  (i.e.   on-peak)  periods.   However, 
determining  precisely  how  fast  to  produce  each  product  during  each  period  of 


each  type  of  day  in  order  to  take  maximum  advantage  of  contractual 
conditions  is  by  no  means  easy. 

C.  Shut-down  Operation 

Given  short-  and  long-term  fluctuations  in  demand,  ArON 
occasionally  has  excess  production  capacity  in  some  regions.   Gaseous 
products  cannot  be  inventoried,  and  inventory  capacity  for  liquids  is  limited. 
Therefore,  it  is  often  necessary  to  put  a  site  into  standby  mode  for  some  part  of 
the  month. 

It  is  here,  in  the  representation  of  site  shut-down,  that  the  principal  use 
of  MIP  arises  in  PRAMSYS.   A  pure  Linear  Programming  (LP)  model  would 
choose  to  shut  a  plant  down  only  during  on-peak  hours,  when  energy  and 
power  are  both  most  expensive.   In  practice,  such  a  solution  would  be 
impractical  for  operational  reasons  (repeatedly  stopping  and  starting 
production  places  unacceptable  stresses  on  equipment  and  requires  an 
excessive  amount  of  operator  intervention.)  While  the  purpose  of  PRAMSYS 
was  not  to  schedule  site  production  day  by  day  or  hour  by  hour,  it  was  vital 
that  the  solutions  be  operationally  feasible.   It  was  therefore  necessary  to 
impose  a  kind  of  loose  parity  between  the  length  of  time  the  site  would  be 
shut  down  during  on,  off,  and  mid-peak. 

D.  Slates:   Discretizing  the  Decision  Space 

Both  energy  and  power  costs  can  be  very  significant.  Since  both  are 
directly  related  to  KW  demand  it  was  clearly  important  to  represent  these 
relationships  accurately.  One  approach  might  have  been  to  use  quadratic 
programming,  but  this  was  rejected  for  several  reasons.   First,  there  are  no 
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commercial  grade  QP  codes  capable  of  handling  MIP  constructs.  Also,  we  had 
at  best  only  an  empirically  derived  quadratic  KW  function. 

We  chose  instead  to  discretize  the  production  rate  space  for  each  site 
into  a  set  of  production  slates.   Each  slate  is  a  vector  containing  a  production 
rate  for  each  product,  and  the  minimum  KW  draw  associated  with  operating 
the  site  to  produce  at  those  rates.   Slates  in  the  set  cover  a  "grid"  of  production 
rates  for  all  products,  closely  approximating  the  actual  KW  production  rate 
function.   The  basic  decision  of  the  model,  therefore,  is  to  determine  how 
long  to  operate  each  potential  slate. 

The  slates  define  the  convex  KW  surface.   This  itself  would  have 
presented  little  problem,  since  LP  is  quite  able  to  represent  a  convex  cost 
function.   However,  the  cost  of  power  (so  called  demand  charge)  is  based  on 
the  maximum  instantaneous  power  demand  over  the  entire  period.    Thus, 
this  cost  would  be  incurred  only  by  the  set  of  production  rates  used  during  the 
period  that  resulted  in  the  greatest  power  draw,  regardless  of  how  long  that 
slate  was  operated.   Representing  this  "if  and  only  if"  relationship  between 
cost  and  production  activity  also  required  the  use  of  MIP  techniques. 

V.    MODEL  FORMULATION 

We  present  here  the  original  MIP  formulation  upon  which  PRAMSYS 
was  based.   Experience  with  MIP  models  drawn  from  this  formulation,  both 
prior  to  and  after  their  application  to  actual  planning  problems,  led  to  a 
number  of  modifications  and  simplifications.   These  are  discussed  briefly  at 
the  conclusion  of  this  section.   In  the  following  section,  we  discuss  our 
approach  for  implementing  the  system  based  on  these  models,  and  experience 
with  the  system. 
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Indices 

i:       1  to  I  index  for  plants 

j:       0  to  J  index  for  slates  at  each  plant  (slate  0  is  plant  shut-down) 

k:      1  to  K  index  for  products 

m:     1  to  M  index  for  customers 

Parameters 

Pjj        =     power  draw  for  jth  slate  at  plant  i  (KW) 

e  i         =     electric  energy  charge  at  plant  i  ($  per  KWH) 

E  J        =     electric  power  demand  charge  at  plant  i  ($  per  KW) 

Cj^jj,     =     cost  of  transporting  one  unit  of  product  k  from  plant  i  to 
customer  m  ($  per  cubic  foot) 

a  jjjj      =     instantaneous  production  rate  of  product  k  by  jth  slate  at  plant  i 
(cubic  feet  per  hour) 

d  i^jj,     =     demand  for  product  k  by  customer  m  (cubic  feet) 

R         =     minimum  run  time  for  any  slate  at  any  plant  (hours) 
T         =     length  of  planning  horizon  (hours) 

Variables 

tjj         =     length  of  time  plant  i  uses  jth  slate  (hours) 

Wj      =     maximal  power  demand  at  plant  i  (KW) 

|l  if  jth  slate  at  plant  i  is  used  at  a  positive  level 
'J  \0  otherwise 

y  ij^jj,    =     quantity  of  product  k  shipped  from  plant  i  to  customer  m  (cubic 
feet) 
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The  objective  function  (1)  in  this  model  is  the  sum  of  energy  costs, 
energy  power  demand  costs,  and  distribution  costs.  Note  that  energy  and 
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power  costs  differ  from  plant  to  plant.  This  is  because  the  contracts  with 

electric  utilities  vary  by  location,  and  furthermore,  each  plant  has  its  unique 

design  and  operating  characteristics.   Note  also  that  the  slates  available  for  use 

at  each  plant,  and  their  costs,  are  uniquely  associated  with  that  plant.   We 

have  chosen  the  fixed  number  J  of  trial  slates  for  each  plant  simply  for 

expositional  convenience. 

The  constraints  (2)  state  that  the  total  quantity  shipped  from  each  plant 

cannot  exceed  the  total  production.   In  practice,  the  inequality  was  extended  to 

account  for  small  quantities  of  beginning  and  allowable  ending  inventories. 

The  constraints  (3)  state  that  the  entire  planning  horizon  is  consumed  at  each 

plant  by  production  time  and  down  time  (recall  that  slate  0  is  the  plant  shut- 
down slate.)  The  constraints  (4a)  and  (4b)  state  that  the  time  t-  that  the  jth 

slate  is  used  at  plant  i,  if  it  is  used  at  all,  must  lie  between  the  conditional 

minimum  R  and  the  maximal  allowable  time  T.    The  upper  bounding 

constraint  in  (4b)  is  redundant  in  the  light  of  constraint  (3);  we  have  included 

it  for  expositional  purposes.   Constraint  (4c)  ensures  that  the  power  demand 
Wj  upon  which  the  power  charge  is  based  equals  the  maximum  of  the  power 

demand  draws  among  slates  selected  by  the  model  for  plant  i.   The  constraints 
(5)  state  that  demand  must  be  met  by  shipments  from  the  plants.   We  note 
that  most  customers  demand  only  one  product.   Thus,  the  total  number  of 
constraints  (6)  is  far  fewer  than  KM. 

The  specific  models  generated  by  PRAMSYS  turned  out  to  be  more 
complex  than  (PAM)  for  several  reasons.    First,  the  model  was  extended  to 
distinguish  among  peak,  mid-peak  and  off-peak  operations  when  the 
electricity  rates  vary  significantly.    Plant  shut-downs  were  modeled  more 
extensively  to  ensure  that  shut-down  periods  occur  contiguously.    Moreover, 
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contracts  with  the  electric  utility  may  be  more  complicated,  involving,  for 
example,  terms  relating  to  differences  in  power  draws  between  peak  and  off- 
peak  periods.  These  complications  were  modeled  by  straightforward 
extensions  of  the  modeling  techniques  used  above.    Finally,  for  complex 
manufacturing  sites  involving  several  interconnected  plants,  the  models 
were  extended  so  that  they  would  choose  the  plant  configurations  as  well  as 
the  slates  for  each  plant. 

Even  without  these  extensions,  (PAM)  is  a  large  scale  MIP  model  of  the 

fixed  charge  variety.   In  particular,  the  power  demand  charges  associated  with 
the  Wj  behave  in  a  manner  similar  to  fixed  charges.   Tricks  involving  cutting 

planes  on  the  plant  objective  functions  derived  from  an  optimal  LP  solution 

proved  relatively  effective  in  causing  the  models  to  produce  good  solutions 
quickly.   A  uniform  reduction  in  size  of  the  demand  charges  Ej  relative  to  the 

energy  charges  ej  also  caused  the  branch  and  bound  to  work  more  efficiently. 
A  second  pass  through  the  MIP  optimization  with  the  best  solution  from  this 
heuristic  as  incumbent  required  far  less  CPU  time  than  that  required  from  a 
cold  start  without  an  incumbent. 

Feedback  from  users  at  the  plants  led  to  an  important  simplification 
that  allowed  the  models  to  be  still  more  rapidly  optimized.    For  the  purposes 
of  monthly  planning,  the  people  running  the  plant  prefer  to  employ  one  slate 
for  each  contract  period  (peak,  mid-peak,  off-peak.)  The  slate  suggested  from 
an  optimal  solution  to  (PAM)  for  each  contract  period  is  the  convex 
combination  of  the  slates  where  the  weights  are  the  fractions  of  the  time  that 
a  slate  is  used.   Since  the  surface  of  the  cost  vs.  slate  function  for  the  plants 
studied  thus  far  has  empirically  proven  to  be  convex,  we  have  been  able  to 
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relax  the  corresponding  MIP  constructs  in  optimizing  the  model.    However, 
MIP  constructs  are  still  required  to  properly  model  shut-downs. 

VI.    IMPLEMENTATION  AND  RESULTS 

PRAMSYS  was  implemented  for  an  IBM  mainframe  computer  using 
the  LOGS  model  generation  language  (see  Brown  et  al  (1986))  and  the  IBM 
optimization  package  MIP/370. 

It  is  important  to  emphasize  that  the  LOGS  model  generation  in 
PRAMSYS  produces  a  family  of  models.   The  precise  formulation  of  a  model 
for  a  specific  region  consisting  of  several  plants  depends  on  the  data  passed  to 
it.   For  example,  depending  upon  whether  a  certain  contractual  element  is 
present  in  the  data,  certain  structures  may  or  may  not  be  present  in  the 
model.   We  reiterate  that  the  model  (PAM)  discussed  in  the  previous  section 
was  merely  the  point  of  departure  for  our  implementation  work,  and  the 
creation  of  a  generator  for  a  family  of  models. 

The  MIP  models  generated  thus  far  for  the  Mid-Atlantic  Region  have 
tended  to  be  quite  large.  As  many  as  1000  slates  for  each  of  several  plants  were 
generated  by  the  Site  Process  Optimization  Protocol  and  included  in  the 
PRAMSYS  models.    Moreover,  the  models  incorporate  upward  of  1000 
customers  demands  over  a  typical  monthly  planning  horizon.    Automatic 
customer  aggregation  procedures  were  implemented,  but  have  not  yet  been 
extensively  used.   The  resulting  models  have  a  few  thousand  rows  and  as 
many  as  10,000  columns.    Using  the  simplifications  and  approximations 
outlined  above,  the  models  are  usually  optimized,  at  least  to  a  close  first 
approximation,  within  a  few  CPU  minutes  on  an  IBM  3083  computer. 
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We  believe  that  the  use  of  PRAMSYS  in  the  Mid-Atlantic  Region  has 
lead  to  shifts  in  the  prevailing  production  and  distribution  patterns. 
However,  as  is  often  the  case  in  real-world  applications,  it  is  difficult  to 
substantiate  this  belief  with  experimental  results,  for  the  simple  and  obvious 
reason  that  PRAMSYS  is  not  run  in  an  experimental  context.    Customer 
demands  fluctuate  from  month  to  month,  and  there  is  no  "control"  process 
to  show  what  would  have  been  done  in  the  absence  of  a  model. 

A  "base  case"  was  run  early  in  the  project,  in  which  PRAMSYS  was 
used  to  second  guess  a  recent  month's  allocation  decisions.   The  model 
solutions  showed  an  increase  in  distribution  costs,  with  a  decrease  in 
production  costs  that  more  than  compensates  for  this  increase.   Overall,  the 
estimate  is  that  PRAMSYS  produces  monthly  production/distribution 
strategies  that  are  several  percentage  points  lower  in  total  cost  than  solutions 
that  would  have  been  obtained  without  it. 

VII.    CONCLUSIONS  AND  FUTURE  RESEARCH 

PRAMSYS  has  proven  itself  to  be  a  useful  and  important  planning  tool 
at  ArON.   Its  success  demonstrates  once  again  that  computer  technology  has 
at  last  reached  a  level  of  development  permitting  mathematical 
programming  models  to  be  implemented  and  effectively  applied  to  business 
planning  problems.  The  success  of  this  project  was  also  due  to  a  felicitous 
blending  of  scientific  skills  and  experience  in  chemical  engineering, 
mathematical  programming,  and  computer  systems  design  and 
programming.    Finally,  the  backing  of  top  management  in  supporting  a 
radically  new  approach  to  planning  was  crucial  to  the  project's  success. 
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PRAMSYS  is  currently  being  extended  for  use  in  other  ArON  national 
regions.   In  this  regard,  experimentation  with  the  Site  Process  Optimization 
Protocol  is  required  for  those  sites  consisting  of  several  production  plants  that 
can  be  linked  in  different  ways.   An  MIP  model  for  calculating  slates  for  these 
more  complex  sites  has  been  developed.   A  related  area  of  future 
experimentation  is  to  link  the  Site  Process  Optimization  Protocol  more 
directly  to  the  PRAMSYS   models  via  price  directed  decomposition  methods 
(see  Shapiro  (1979).  The  idea  would  be  to  occasionally  use  shadow  prices  from 
the  mathematical  programming  model  to  price  out  slates  produced  by  SIPOP, 
and  select  new  slates  for  the  PRAMSYS  model. 

Once  models  for  all  regions  have  been  developed,  the  intention  is  to 
construct  a  longer  range,  national  model  for  strategic  planning  purposes.   The 
types  of  problems  to  be  addressed  by  such  a  model  include  contract 
negotiations  with  customers  and  electric  utilities,  long  term  plant  shut- 
downs, and  economic  evaluations  of  new  markets. 

Moving  in  the  other  direction  with  respect  to  time  and  scope,  a  new 
project  is  underway  to  convert  the  production  planning  sub-model  in 
PRAMSYS  to  a  production  scheduling  model.   The  reader  may  have  noted 
that  the  model  (PAM)  selects  an  optimal  combination  of  slates,  but  makes  no 
attempt  to  schedule  the  sequence  in  which  they  should  be  used.    In  the  short- 
term  when  one  considers  distinct  production  periods  with  varying  demands 
on  the  plant,  and  recognize  that  inventory  storage  for  gas  is  extremely 
limited,  the  sequencing  of  slates  becomes  important. 

Finally,  generalizations  of  the  models  developed  for  PRAMSYS  should 
be  applicable  to  other  process  manufacturing  industries.   The  underlying 
principle  in  perfornMng  modeling  research  in  this  area  is  to  better  understand 
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how  to  imbed  process  control  optimization  models,  which  provide  an 
instantaneous  prescription  for  the  plant,  in  one  or  more  mathematical 
programming  models  for  planning  and  scheduling.    A  central 
methodological  problem  is  how  to  pass  from  the  essentially  "instantaneous" 
world  of  process  control  to  the  world  of  operational  scheduling  and  control, 
in  which  sequences  of  discrete,  on /off  events  associated  with  changeovers 
and  setups  are  crucial.    We  believe  the  models  in  PRAMSYS  are  an  important 
step  in  this  research  direction. 
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